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The scattered field produced by a plane wave incident on a perfectly conducting sphere 
coated with a thin layer of material w-ith large complex index of refraction is considered. 
It is shown that, for certain relations involving the thickness of the layer, and its permittivity 
and permeability, the scattering problem is equivalent to the problem wherein the total 
electric and magnetic field components must satisfy a general impedance boundary condition 
on the outer surface. With this simplification the backscattered field is obtained in terms 
of the geometric optics portion and the diffracted field portion for small wavelengths. 

1, Introduction 

Previous work on scattering by spheres comprised of concentric layers of various materials 
has been carried out by Aden and Kerker [1951], and Sluirfman [1954], who obtained the exact 
Mie series for the various cases they considered. However, they did not perform an asymptotic 
evaluation of the fields for the small wavelength case. 

In this paper the investigation for small wavelengths is considered foi; the case w^here the 
particular diffracting body is a perfectly conducting sphere coated with a lliin layer of large 
complex refractive index. First the exact Mie series is given for the scattei'ed field. Using 
asymptotic methods it is shown that the effect of the coating can be replaced by an equivalent 
impedance boundary on the outer surface. In so doing, certain r(»strictions ai*e placed on the 
thickness of the coating and the argument of the index of refraction. 

A similar problem, namely that of a radial electric dipole, outside a concentrically' strati- 
fied spherical conductor was treated by Wait [1956]. An approximation obtaining the imped- 
ance boundary condition was also discussed by Wait. 

2. Mie Series 

Consider a perfectly conducting sphere with radius b covered by a homogeneous layer with 
uniform thickness 8 of a material characterized by its electric permittivity e and relative mag- 
netic permeability /z' — which may be complex- — on which is incident the plane electromagnetic 
wave 

where 

, 27r oj . .^,2 

A C 

r;o=(eo/Mo)^''^, 

and €o and fiQ are the permittivity and permeability of the surrounding medium, taken to be 
free space. The wave number of the coating is given b}' ki and is related to the index of refrac- 
tion by the relation 

h=Nk. 



1 The work described in tliis paper was carried out for the Air Force Ballistic Missile Division under subcontract to the Conductron Corporation 
as a part of the progiani being carried out linder Chrysler Corporation Prime Contract Ar04(694)-25. 
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The general Mie series solution for a sphere with a concentric shell is given b}^ Aden and 
Kerker [1951] and the case where the inner sphere is a perfect conductor by Sharfman [1954]. 
Using these results the far scattered field is thus found to be, in spherical components form, 

Er=Hr=0 

^ yl/n{ka)An—ypn(ka) 

7 ^ i/n{ka)Bn—yly'n{ka) .^. 

^^ ^'^\ka)B,-i',''\ka) ^ ^ 

and \l/{ka) and ^(ka) are related to the Bessel functions b}^ 

with j/ = /i+K- 

The radius of the outer sphere is given by a = b^8. The functions A^. and Bn are given by 

Wlttl 

a= f ^^^ (ha) f i^> (A:i6) - f ;?^ {h b) f i^> (tia) . (5) 

3. Impedance Boundary Condition Approximation 

It will be shown that for large ka, N, kb, and certain other conditions, the exact Mie series 
for the scattered field can be approximated by a Mie series for which the total field on the sur- 
face of the sphere satisfies the generalized impedance boundary condition, Leontovich [appen- 
dix to the Logan and Blacksmith collection of Fock papers, 1957] 

E-(n.E)n-r777o-^nXH (6) 

where the relative impedance r? is a function of ^\ e, and 6. 

The procedure is centered around the asymptotic approximations to the coefficients A^ 
and Bn, First the coefficient Gn is considered. It is written in the following two forms 

C,='^^h^M {2j,{ha)m'\hb)-2j,{k,b)m'\k,a)} 



=^ ^ha kib {Hi'' {ha)m'\hb)-m'\hb)m'\kia) } (7) 

where z^=7i+/2- Before considering the asymptotic evaluation of the Hankel functions for 
large order, the following restriction will be placed upon the argument of N, 

0<argAr<^-^0' (8) 
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where the positive angles unci 0' satisf3^ the inequalities 

0»|iVH|-2/^ 

0'»\Nkb\-K (9) 

The lowTT bound on the arguineMt of iV insures that the following inequalities hold, 

\ha-v\>v'^' 

\hh-v\yv"^ 

for the complete range of n of interest. Thus the Debye asymptotic expansions for the Bessel 
and Hankel functions may be used. Hence from Watson [ch. 8, 1958], the coefficient C„ has 
the following asymptotic form for large Nka and Nkh, 

2 

V( — '^ smh7i)(— 1 smh72) 
where 

v=Nka cosh 7 
v=Nkb cosh 72 
with 

7i=q;i + ?:/3i, and i2=0L2^-i^2- 



;} "" 



Both ai and a^ are negative and both /3i and ^2 lie between and 7r/2. The functions i^ 
are given by 



i = j^(tanh7i-7i)'y^ .^2) 

2 = ^(tanh 70— 72) J 



Expression (10) is the leading term of an asymptotic expansion, provided that condition (8) 
holds. If the argument of A^ is the order of, or less than 0, then the Airy integral approx- 
imation would liave to be used for the Bessel functions when v/{Nka) or v/{Nkb)^l. If the 
argument of N is the order of, or greater than 7r/2 — 0', then expression (10) would have to 
be modified by the addition of extra terms when p/(Nka) or v/(Nkb) is in the vicinity of ~1.5i. 

Without going into their derivation, certain properties of the functions associated with the 
right-hand side of eq (10) will be given. First, it can be shown that the real part of 1/^2 — ^Ai is a 
positive function and jnonotonically increases with increasing real v. The imaginary part of 
\^2— 'Ai is negative and monotonically increases with increasing real v. The real part of sinh 71 
and sinh 72 are both negative, whereas their imaginary parts are both positive. 

From (3) and (10) the asymptotic behavior of An for large A^, ka, kb is given by 



and similarly for B^ 



N 
An 7 coth (lAo-^i) sinh 7i (13) 



B„ ^ tanh {^2- ^1) sinh 71. (14) 



111 particular for n<^<! \Nkb\, these become 

A~-f coth (^.-w{l + 0(^y+0(^,)} (15) 

5„^-^tanh(^.-^,){l+0(^J+0(^,)} (16) 

where the term wliich is the order of d/{kia^) arises from the consideration of Iiigher order terms 
in the asymptotic expansion for A^ and B^. Finally one obtains 

615 



^-*.=-*.'+*'S(^)+»[K5)'} <'« 

Hence for the range of n, 0<Cn<^kaM where M is a number unspecified as yet but lying in the 
range 1<CM<C<C |A7^|, the following approximation holds, 

V'2-'Ai ih8 (18) 

provided that d is sufRciently small such that 

kx8 



M' 



2N' 



«L (19) 



Thus for the range 0^n<^Mka, the constants An and Bn are effectively independent of n. 
Hence set 

Bn=-ir}. (20) 

Using the relation An'^ — (Bn)~^ the set of coefficients a^ and bn of the Mie series (eq (21)) can 
be approximated b}^ the following, 



\pn{ka)—iv\pn{l^a) 
^i''{ka)-iv^'n''\ka) 






ir}\l/n{ka) + \pn(ka) ,^2) 



ivtn'\ka) + ^'n'''{ka) 



But these are the corresponding coefficients that one would obtain for a boundary condition 
of the type given by eq (6) [see the appendix of Hiatt et al., I960]. 

The remaining problem is to consider the behavior of the coefficients for n^Mka. It will 
be assumed that M is sufficiently larger than unity so that the Debye asymptotic expansions 
may be used for \pn(ka) and ^l^^(ka) for ny>Mka, Explicity M must be such that M— 1>> 
(ka)~^^^. Thus one obtains 



exp {j^(tanh7— 7) 



Uka)^ '''''"-'^:^ '" (23) 

Vsmh 7 

^<,^^^^^expi-^;ganM:^ (24) 

Vsinh 7 

when v=n-\-%=ka cosh 7 and 7 is positive and real. 

Using the above approximations the coefficients a^ and bn become 

r isinh7i + (jLt7iV) sinh7 tanh (1^2— 1A1) n o r^- i. \ /ok\ 

, r i sinh Ti + (iV/M') sinh 7 coth (^^2-^1) "| „„„ p ,., . ,oa\ 

For n sufficiently greater than ka^ the exponents in the above expression are large and nega- 
tive. Hence provided that the denominators of the expressions given by (25) and (26) do not 
become too small, the coefficients a^ and bn for n sufficiently greater than ka, are negligible. 
In consideration of the denominators it can be shown that since the real part of (1^2—^1) is 
positive and gi*eater than the real part of —iki8, then the real part of coth (1^2—^1) is positive 
and with respect to magnitude, is greater than the imaginary part of coth (1^2— ^Ai)- In particular 
as n approaches 00^ coth (^2—^1) approaches unity. A similar anatysis holds for tanh (^2— ^i)- 
Hence if 

|arg(iV/M')l<j (27) 
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then 

Real {(//N) taiih (^l^2-^l)}>0 

Real [{N/fi') coth (^2-\^i)}>0. 

In addition since the imaginary part of sinh 71 is positive^ then 

Real {—i sinh 71 } >0. 

Hence it is easily seen that the real parts of the denominators of the rigiit-hand sides of expres- 
sions (25) and (26) are greater than zero. Condition (27) can be relaxed to permit a greater 
range of permissible values for the argument of N/^i^, but only after considerable and extensive 
analysis. 

Considering the corresponding problem for the coefficients of the form given by (21) and 
(22) one obtains for ny-Mka 



bn^l —T- r— j ) exp 2v (tanh 7—7). 

\ iry— smh 7 / 



Provided that 77 is bounded away from the imaginary axis, these behave like expou(Mitials with 
large negative exponents. 

Thus the Mie series for the coated spliere may be approximated by the Mie series for the 
boundary condition on the outer surface given by eq (6) with 

r;=~'^tan(^i5). (28) 

4. High-Frequency Backscattered Field 

For the evaluation of the bis ta tic and forward-scattering far fields in the high-frequency 
range, it is best to express the scattered field in terms of the Debye potentials. However, 
for the case of backscattering one can work directly with the expression for the electric intensity. 

From eq (1) the backscattered field is given by 

E'=Eo~:(-i) g (-l)"(«+l/2)(a„-W. (29) 

Simplification is achieved for the sphere with a thin coating of large refractive index if the 
approximations to the coefficients a^ and bn given by eqs (21) and (22) are used. 

Combining the two coefficients and using the Wronskian relations for the Bessel functions, 
one obtains 

an—hn==^ ~\ Ln + i/2{'n)Ln-i-l/2{^/v) (30) 

where 

Treating the summation over ti as a residue series, the summation is replaced by a contour 
integral C taken in the clockwise direction around the poles at 1^=1/2, 3/2, . . ., giving 

kr V.2' 2rj JcCOsvTrL,{'n)L,{l/v) J 

The contour integral can be deformed using the Watson transform technique to give a line 
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integral extending from the fourth quadrant through the origin to the second quadrant in 
the V plane, plus a residue summation produced by the poles of Lp{r)) and Lr,{l/r}) in the first 
quadrant. The details of this well-known technique will not be discussed here. 

An asymptotic evaluation of the line integral can be performed in the same manner as the 
technique of Scott [1949]. The combination of the asymptotic expression for the line integral 
and the term involving (ao—bo) in the expression (32) gives the geometric optics contribution 
to the backscattered field, namely 



^-' g.o. 



^(hr^{:i)V-Wa^ ■ ■ y'-''--^^- (^^> 



The residue series corresponding to the ^ ^creeping waves'' or diffracted field contribution to the 
backscattered field is 

(34) 

where vi and Vs ^^i'^ the zeros of Lp(rj) and Luil/rj) respectively. Since the dominant contribution 
arises from the zeros which are close to ka, the following asymptotic approximation is made, 

^i%2ika)^-im'^'wi{t) (35) 

^:^-'h2(ka)^im-'''w[(t) (36) 



where 

m 



=(— j and j'^te+m^. 



Here the notation of Fock [1946] is used for the Airy integral Wi(t). It is related to the notation 
of Miller [1946], which is more useful for numerical results, by the following expression, 

w,{t)=^[BS)+iA,{t)]. 

For further reference on the subject, one should consult Logan [1959]. Excluding the case 
where r;=0, which corresponds to the perfect conducting sphere, whose results are well-known, 
the residue series can be approximated by 

^0- Is -^^ [(^'+—) (Mti)yT'-f:-~'^^[(\+—^ (wdtswYX (37) 

m \^f^i cos VfT LV my ' ' J l^cos VsTt [_\r my J J 



kr 
where ti and ts are respectiveh^ the roots of 

w[{t) 



--irim (38) 



^i(^) 
and 

^=rW,. (39) 

Wi(t) 

The values of w[ {t)jwi (t) have been plotted by Logan and Yee for values of t in the first 
quadrant. From this the first few roots of the eqs (38) and (39) can be obtained. 

The total backscattered field is obtained by combining expressions (33) and (37). 

The reduction of the original problem (that of scattering by a perfectly conducting sphere 
with a thin coating of material of large index of refraction) to the equivalent problem (where 
the boundary condition is expressed in terms of an impedance boundary condition) greatly 
reduces the effort in obtaining the small wavelength asymptotic expression for the far scattered 
field. A similar technique should be applicable to other coated bodies with large radii of 
curvature. 
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